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1. INTRODUCTION 


Adsorption of metals in soils, sediments and aquifers is a 
complex phenomenon. Metal ions are charged species which may form 
various soluble complexes both with inorganic and organic ligands. 
The formation constants for the inorganic ligand-metal complexes are 
well known for many complexes and can be used to calculate the 
distribution of the metal over its various forms in solution. However, 
for complexes between metals and natural dissolved organic ligands, 
like fulvic acids, it is not possible to calculate the speciation. This is 
due to the absence of reliable physical chemical models that can 
account for metal binding to natural organic matter. This seriously 
hinders the advance of research related to metal speciation in the 
environment. 

The presence of soluble metal complexes often reduces the metal 
adsorption compared with the absence of these dissolved complexes. 
This is the case when the affinity of the metal complexes for reactive 
sites on the solid matrix is much less than that of the free metal ion. 
Examples are metal complexes with ligands like chloride or citrate [1]. 
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Another important factor influencing metal adsorption to natural 
sorbents is the competition between a metal ion of interest and other 
species for the same adsorption sites. An example is the influence of 
zinc and calcium on the cadmium adsorption to soils [1]. In practical 
pollution situations a mixture of metals and other ions is often involved. 
Competition may, in such cases, notably reduce the adsorption of a 
weakly adsorbed metal ion like cadmium, causing enhanced mobility 
and bioavailability [2]. The retardation of a metal in groundwater is 
thus strongly a factor of the solution phase composition, which will 
vary with respect to, for example, the position in the pollutant plume. 
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FIGURE 1. Calculated Adsorption Isotherms of Cadmium on Iron 
Oxides for Various CI Levels, Represented with a Total 
Cadmium Solution Concentration (a) and a Cd ** Activity 
Scale (b). The cadmium adsorption constant (log K = 6.8) 
has been taken from Van Riemsdijk et al. [24]. The Stern 
layer capacitance is set at 2.23 and only one protonation 
reaction is assumed (log K,, = 10.7). 
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In Figure 1, the effect of chloride on cadmium adsorption is 
illustrated. When the adsorption is plotted as a function of total metal 
concentration in solution a whole range of curves results, the actual 
shape depending on the chloride concentration in the solution. These 
curves coalesce approximately into one curve when the adsorption is 
plotted as a function of free metal ion activity in solution, (Figure 1b), 
provided that the adsorption of the chloride complexes is negligible and 
that "competing" species are present at constant solution activity. The 
behavior shown in Figure | has been observed for cadmium adsorption 
for a whole range of different soils [1]. 

It is well known that pH strongly influences the adsorption of 
metal ions in natural systems. This effect is caused firstly by a change 
of the surface potential of the reactive solid phases and secondly by a 
change in the competition between protons and metal ions for the same 
adsorption sites. A change in the surface potential affects metal 
adsorption through a change in electrostatic repulsion or attraction. 

Apart from changes in solution composition, changes in solid 
phase composition will also strongly influence the interaction of metals 
with the solid matrix. The question arises how to deal with the above 
described complexity. A ‘practical’ approach would be to measure the 
interaction with a natural solid matrix of interest for a range of solution 
conditions that are considered relevant for the situation that is to be 
dealt with. One can then try to model the obtained result with a chosen 
empirical model using curve fitting techniques to obtain the model 
parameters. Although this approach may give a good description for 
the situation studied, nothing can be said with any confidence for a 
situation that is outside the range of conditions that have been used in 
the experiments. Nor can the information obtained be transferred to a 
solid matrix with a different composition of reactive phases. For the 
purpose of environmental policy, soil chemical knowledge is often 
required that is of more general validity. 
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2. CRITICAL EVALUATION OF EXISTENT 
(CHEMICAL) METHODOLOGIES DEALING WITH 
THE COMPLEXITY OF SOIL SYSTEMS 


The above indicates that natural variation in both the solution 
composition and the composition of the solid matrix may lead to a very 
large variation in the distribution behavior of one and the same 
chemical for different locations. It is thus far from simple to obtain 
information that has general validity. Several approaches can be found 
in the soil chemical literature that are attempts to obtain this type of 
information. 

One approach is to discriminate the total content of a chemical 
over various operationally defined fractions in the solid phase. This 
solid phase "speciation" is obtained by employing sequential extraction 
procedures. For metal speciation in sediments, the following fractions 
are "defined" by sequential extraction [3-7]: 


1) exchangeable metals plus metals associated with calcium 
carbonate, 
2) metals present in easily reduced phases, mainly metals 


associated with manganese oxides, and partly also associated 
with amorphous iron oxides, 


3) metals associated with amorphous and poorly crystalline iron 
oxides, 

4) metals associated with organic matter or present as metal 
sulphides, and 

5) a residual fraction. 


Sequential extraction with different extraction agents by 
definition leads to various metal fractions. Whether the fractions 
obtained indeed represent the various pools is questionable. Gruebel et 
al. [8] tested a sequential extraction procedure for arsenic and selenium 
using well characterized minerals. Their work shows that the various 
phases may not always be identified correctly. 

Identifying operationally defined fractions in itself does not lead 
to the possibility of predicting the interaction between a chemical and 
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various soils, sediments and aquifer materials. One step in that 
direction is the work of Unger and Allen [7] who combined the 
measurement of the partitioning of a metal over the solution phase and 
the sediment with a sequential extraction procedure. In this way, they 
were able to "speciate" the amount (ad)sorbed as a function of the 
loading of the sediment with the metal. They estimated a binding 
capacity and an average equilibrium "constant" for each of the five 
fractions that were discriminated in this study. How this binding 
capacity is related to soil or sediment type was not a subject of their 
study. 

The relation between soil characteristics and the observed overall 
interaction behavior of a chemical can be studied by measuring the 
partitioning of a chemical over solution and solid phase for a whole 
range of different soils and sediments [9,10]. Although the 
methodology in itself is promising, there are a few serious problems 
that may explain why the results published so far are often 
disappointing. The first problem is that the soil characteristics 
measured are poorly (operationally) defined and are not necessarily 
good estimators for the various types of reactive surfaces present. 

The second, very important, problem is that the chemical 
composition of the aqueous phase from which the adsorption is studied 
generally varies for the different soils. Although the same electrolyte 
may be added to all soils, the resulting aqueous medium is not the 
same. An obvious and very important factor is the pH of the aqueous 
phase, which is buffered by the soil particles, resulting in different pH 
values of the electrolyte solution for different soils. Because the pH 
strongly influences most chemical interactions in soil, the result of such 
a study is sometimes that the pH is the soil characteristic responsible 
for differences in behavior between various soils [10]. This suggests 
that different soils, when brought to the same pH, should show the 
same adsorption behavior. This is obviously not true. Another example 
of unintentional differences in the resulting composition of the aqueous 
phase is the presence of a variable amount of competitor species. 
Studying cadmium adsorption for a range of soils will, in general, lead 
to a release of a variable amount of zinc ions. The resulting cadmium 
‘adsorption i is not only a function of the various reactive surfaces that 
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are present and the resulting pH, but also depends on the degree of 
competition with zinc ions [11]. 

The third problem is that soils as sampled may already contain 
a varying amount of the adsorbed species of interest. To be able to 
make a good comparison between the different soils, the desorption 
behavior of the species of interest should be studied. Studying 
desorption is often extremely difficult, especially in the case of a high 
affinity adsorption, i.e. a steep ad(de)sorption isotherm. Recently a 
promising new methodology for studying phosphate desorption has been 
developed [12]. 

Another approach in trying to obtain results that are of a more 
general validity for a wide range of soils is to study the interaction of 
a chemical with various individual soil components as a means to 
interpret the difference in behavior between various soils or sediments. 
An example of this approach is the work of McLaren et al. [9] who 
studied adsorption of cobalt with different soil components and also 
with different soils. At present, this procedure is also not very 
successful mainly because for two reasons. 

The first reason is that the fundamental knowledge about the 
physicochemical interaction mechanisms between species and the 
various soil components is, at present, not well enough established. 
McLaren et al. [9] used empirical equations to describe the measured 
interaction with soil components. The measurements with the soil 
components were done for a very limited range of solution conditions, 
at one pH only, and the empirical model used does not allow for 
estimating the interaction behavior for solution conditions other than 
those of the measurements. The adsorption measurements with the 
different soils were done with one electrolyte solution without control 
of pH leading to different pH values for the different soils studied. 
Comparison between these adsorption measurements and those of soil 
components based on one pH measurement is bound to fail. The 
second important problem with this methodology is the need to 
establish the content and nature of the various reactive surfaces present 
in a soil sample. 

For the development of general applicable soil chemical 
knowledge that can be used to interpret the bio-availablity, toxicity and 
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mobility of pollutants, it is essential to have insight in the fundamental 
physical chemical interaction of species with soil components. Metal 
(hydr)oxides of aluminum, iron and manganese and amorphous 
aluminum silicates, together with natural organic matter, are very 
important reactive surfaces with respect to metal adsorption in soils, 
lakes and aquifers. In the following, we will discuss physical chemical 
adsorption models for metal oxides and organic matter together with a 
discussion on how this information can be transferred to natural 
systems. 


3. ADSORPTION MODELS FOR METAL OXIDES 


Metal (hydr)oxides of iron, aluminum and manganese are 
important reactive surfaces in soils, aquifers and sediments with respect 
to interaction with positively charged species like, H*, Al**, Cd”, Zn”, 
Pb**, and Cu” and with negatively charged species like phosphate, 
arsenate, sulphate, selenite, borate, bicarbonate and fluoride. Organic 
species may also adsorb on metal oxides [13]. The adsorption of 
species on oxides is strongly dependent on the pH due to the variable, 
pH dependent, surface charge and potential of metal (hydr)oxides and 
to a pH dependent speciation of the adsorbate. 


4. RELATIONSHIP BETWEEN SURFACE CHARGE 
(POTENTIAL) AND pH 


Metal oxides are amphoteric surfaces which may be either 
positively or negatively charged depending on the pH of the system. 
The pH at which the metal(hydr)oxide is uncharged is called the point 
of zero charge, pzc. The pzc of metal (hydr)oxides differs considerably, 
ranging from 2 for silica to around 10 for some iron oxides. The 
primary surface charge, in the absence of specific adsorptions of ions 
other than H*, depends on the pH and on the electrolyte concentration 
or ionic strength. The most widely used model to interpret the variable 
charge characteristics is the so called two-pK model [13-18]. In this 


8 Riemsdijk and Hiemstra 


model, it is assumed that there is only one type of reactive group 
present on the surface of a metal (hydr)oxide that may react with 
protons according to the following two reactions, 


S07 +H, + SOH*; Ky (1) 


SOH” +H, = SOH; ; Ky (2) 


where SO", SOH°, SOH,* represent three different surface species 
belonging to one and the same reactive surface oxygen group, and [H,*] 
is the proton concentration near the plane of adsorption which is related 
to the concentration in solution, [H*], via, 


[H,] = [H*] exp(-F y, / RT) (3) 


where y, is the surface potential. In order to be able to model the 
variable charge characteristics, an additional expression is required that 
relates the surface charge and the surface potential. This relationship 
is provided by assuming a certain double layer model. In the "oxide" 
literature, different approaches can be found. Sometimes only a diffuse 
Gouy Chapman double layer is assumed [19], or it is assumed that ratio 
between surface charge and potential is constant, leading to the so 
called constant capacitance model [13,14,17], or a combination of a 
Stern layer with a Gouy Chapman double layer often indicated as the 
Basic Stern model is used [20-25], or a model where three planes near 
the surface plus a diffuse double layer model is used which is called the 
triple layer model [15,18,26]. The basic charging reactions in all these 
approaches are equal. They only differ in the relation(s) used to relate 
the surface charge and the surface potential. A disadvantage of the 
constant capacitance model is that it cannot be used to describe the 
charging at variable salt level because the "constant" capacitance is in 
practice a function of the salt concentration. A disadvantage of the 
triple layer model, as used in literature is that in addition to the 
reactions with protons, pair formation with simple electrolyte ions, such 
as sodium and chloride, has to also be included in the model, and this 
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increases the number of adjustable parameters. These variable charge, 
variable potential models have been successfully used to describe 
surface charge-pH curves as a function of salt level for various metal 
(hydr)oxides. The model parameters, e.g. the two proton affinity 
constants and the Stern layer capacitance are treated as adjustable 
parameters. The proton affinity constants may differ considerably for 
various metal (hydr)oxides. 

An alternative for the two-pK model is the one-pK model [20- 
25,27]. This assumes that the basic charging reaction can be described 
with the following reaction, 

SOH" +H; = SOH” ; K, (4) 

The surface charge is dependent on the ratio between the two different 
species, belonging to the same surface group. It follows from this 
model that the log Ką equals the pzc of the metal (hydr)oxide. The 
advantage of the one-pK model is its extreme simplicity, the fact that 
it has one parameter less than the two-pK model and that the proton 
affinity constant follows directly from the pzc. This one-pK model 
gives as good a description of the available data as the two-pK model 
for most metal oxides, with the exception of silica. 


5. THE MUSIC MODEL 


Both the one and the two-pK modeling approaches presume the 
presence of only one type of reactive group. However, it is well known 
from infrared analysis that different types of surface oxygen groups are 
present on metal (hydr)oxides [28-30]. Surface oxygen may be singly, 
doubly or triply coordinated to an underlying metal ion. It is to be 
expected that the reactivity of these oxygen groups towards proton 
adsorption differs considerably. Recently a new model, the MUSIC 
(MUti SIte Complexation) model, has been developed that takes these 
three different types of surface groups into account [21,22]. Each type 
of surface group may react with two protons. The affinity constants for 
the proton adsorption reactions for the various groups of a range of 
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different metal (hydr)oxides are predicted by the new model, leading 
to the possibility of predicting the variable charge behavior of metal 
oxides. It predicts that the difference between the log K values for two 
consecutive proton adsorption reactions for one type of surface group 
is very large (around 14 log K units). The one and two-pK models are 
special cases of this more general new model. The predictions are in 
good agreement with existing experimental data [22]. Another result 
is that the quite different shape of the charging curve for silica 
compared to other oxides is predicted by the MUSIC model. It also 
follows from the model that different crystal planes of the same metal 
(hydr)oxide may exhibit a quite different charging behavior, due to a 
different composition with respect to the various types of reactive 
surface groups. This may have important implications for the specific 
adsorption of species other than protons on metal hydr(oxides). For 
gibbsite (Al(OH),) it is predicted that the dominant crystal plane (001) 
is uncharged for pH values below ten, whereas a negative surface 
charge develops for pH values higher than ten. In the normal pH range, 
only the edge faces will develop a (positive) surface charge. These 
predictions are also in accordance with experimental information [22]. 
Specific adsorption of anions on gibbsite will thus probably only occur 
on the edge faces. 

The prediction of the proton affinity constants is based on the 
following formula [21]: 


log K = A-~Bn — (5) 


where A and B are constants, n is the number of metal ions that 
coordinate with a surface oxygen group (n = 1, 2 or 3), L is the 
distance between an adsorbed proton and the metal center and the 
formal bond valence v, as introduced by Pauling [31] equals, 
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ps (6) 


where Z is the charge of the metal in the metal (hydr)oxide and CN is 
the coordination number of this metal. The constant A differs whether 
a proton reacts with a ‘naked’ surface oxygen or with a surface 
hydroxyl. The constants A and B are obtained by considering the 
protonation of solution monomers for which log K and v/L are known. 
For these monomers, n is always one. The model is further calibrated 
based on the protonation of a singly coordinated surface hydroxyl as 
present on gibbsite (Al(OH),). The log K value for this reaction can be 
derived directly from experimental data because the charging behavior 
of gibbsite is determined by only this one reactive group over a wide 
pH range [21,22]. It turns out that the affinity for proton adsorption on 
a surface is approximately 4 log K units higher than the corresponding 
reaction for the dissolved monomer. This difference may arise because 
the spatial configuration of the oxygen ligands on a surface is in general 
quite different from the situation for a dissolved monomer. This 
difference in configuration will lead to a different enthalpy and entropy 
of adsorption for the surface compared with the reactions in solution. 
The formal bond valence v can be seen as the positive charge 
of the metal center that is contributed to a ligand. For iron and 
aluminum (v = 3/6 = 1/2). The concept of v also facilitates a proper 
bookkeeping of surface charges. This is illustrated in Figure 2, where 
the concept is applied to three different dissolved iron species. The 
positive charge of the central metal ion has been equally distributed 
Over the surrounding ligands. Note that in this case a positive, neutral 
or negative species can be formed by simply changing the ratio of the 
Coordinated hydroxyl and water molecules. This is the basis of the 
‘One-pK metal oxide model. 
= The proton affinity constants for the three different types of 
‘surface oxygen groups that may be present on the surface of iron and 
‘aluminum (hydr)oxide as predicted by the MUSIC model are given in 
‘able 1. Assuming that the entropy contribution to the protonation 
oa tion is mainly derived from the change in hydration entropy, and 
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Fe(OH) (OH), <> Fe(OH) (OH); <> Fe(OH) (OH), 
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FIGURE 2. Schematic Representation of Three Dissolved Fe-hydroxy 
Complexes. 


applying the Born treatment of hydration energy [32] to the protonation 
of solution monomers, and using the formal bond valence in case of 
surface reactions, the entropy and thus also the enthalpy (because the 
free energy can be derived from the log K,,,,.) of the surface protonation 
reactions can be estimated [32]. The results are given in Table 1. It 
follows from Table 1 that the difference in log K between the first and 
the second proton adsorbing on the same type of oxygen group is 
approximately 14 log K units. The difference between corresponding 
protonation reactions for different types of groups (singly, doubly, triply 
coordinated) is also very large. Note that the protonation for a doubly 
coordinated surface oxygen for iron and aluminum (hydr)oxides is 
equivalent with the formulation of the classical two-pK model. The 
difference being that in the two-pK model the adsorption constants are 
treated as fitting parameters, whereas here they are a priori predicted. 
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The protonation of a singly coordinated surface hydroxy] for these types 
of metal (hydr)oxides is equivalent with the formulation of the classical 
one-pK model. 


TABLE 1. Predicted log K for Surface Protonation Reactions at 298 °K 
[21] and for Related Thermodynamic Quantities [32]. 


Reactant v log K G H S 
kJ mol” kJ mol! J mol’ K' 

AlO -3/2 23.8 -135.7 -106+2 101+6 
AIOH -1/2 10.0 -57.0 -42+2 50+6 
Al,O -1 12.3 -70.1 -47+2 77+6 
Al,OH 0 -1.5 8.6 15+2 22+6 
ALO -1/2 22 -12.5 2+2 50+6 
ALOH +1/2 -11.6 66.1 65+2 -4+6 
FeO -3/2 24.5 -139.7 -116+4 79+4 
FeOH -1/2 10.7 -61.0 -52+4 30413 
Fe,O -1 13:7 -78.1 -61+4 56413 
Fe,OH 0 -0.10 0.6 2+4 4+13 
Fe,O -1/2 4.3 -24.5 -16+4 30413 
Fe,OH +1/2 -9.5 54.2 48+4 -20+13 


The log K for the first protonation step of a singly coordinated 
oxygen, formal charge -1.5, is so high that a fully dissociated singly 
coordinated group will not be present to any significant extent on the 
surface of an iron or aluminum (hydr)oxide in contact with water. The 
log K for the protonation of a triply coordinated hydroxyl group is so 
small that a triply coordinated surface water molecule, formal charge 
+1.5, can also be considered to be nonexistent for these (hydr)oxides. 

The outcome of the MUSIC model implies that different crystal 
planes may exhibit completely different properties. For clays it is 
common practice to consider that different planes may have quite 
different properties. However for metal oxides up till now it has been 
common practice to treat a metal oxide as if it has only one type of 
reactive group on the surface. The difference in properties of different 
faces that are predicted by the MUSIC model are first illustrated for 
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gibbsite. Gibbsite crystals are thin hexagonal platelets. The dominant 
001 plane exhibits only doubly coordinated groups on the surface. It 
is predicted that this plane is essentially uncharged below pH 9.0. 
Above pH 9.0, a negative surface charge will develop on this plane. 
Both singly and doubly coordinated groups occur on the edge faces. 
The singly coordinated groups are on these faces responsible for the 
charging behavior over the whole pH range. The doubly coordinated 
groups hardly contribute to the charging of the edge faces even at high 
pH, because the dissociation of the doubly coordinated surface 
hydroxyls is strongly suppressed by the negative surface potential that 
is caused by the dissociation of the singly coordinated groups. The 
effective charging behavior of a group on a surface is thus not only a 
function of its proton affinity constants, but depends also on the 
presence of other types of reactive groups on the same mineral face. 
The overall particle surface charge above pH 9.0 is determined by 
contributions of both the edge faces and the planar surfaces. These 
predictions are in very good agreement with experimental data [22]. 
These observations also explain why the slope of the charging curve of 
gibbsite below pH 9.0 at a certain salt level is relatively low when the 
charge is expressed per unit area BET surface. 

Another metal oxide for which the surface structure is well 
established is goethite [33] as prepared according to Atkinson’s method 
[34]. Three different crystal planes are present [33]. The total surface 
area, for particles with a specific surface area of around 100 m’/g, is 
mainly determined by the 100 and the 010 face, the contribution of the 
two faces to the total surface area being roughly equal. The 
composition of these two faces with respect to the reactive groups is 
quite different [35,22]. The 010 face exhibits only singly and doubly 
coordinated groups. The outcome of the MUSIC model predicts that 
the charging behavior of this face is almost exclusively determined by 
the singly coordinated groups. The charging of this face can essentially 
be calculated with the classical one-pK model. The pzc of this face is 
predicted to be 10.7. At the 100 face singly, doubly and triply 
coordinated groups are present in equal densities, 3.3 sites/nm* each. 
The singly and triply coordinated groups are the reactive groups. The 
pzc of a hypothetical face with only triply coordinated would be at 
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pH=4.3. The combination of the two reactive groups, where the 
charging of each can be described with a single affinity constant, leads 
to a predicted pzc of 7.5 for the 100 face. Around the pzc of the 100 
face of goethite, the singly and triply coordinated groups are nearly 
fully but oppositely charged. This may strongly influence specific 
adsorption of ions at this face. 


Goethite 0.1 M 


Š 


Surface charge (mC/n?) 


FIGURE 3. Calculated 6,(pH) Curves for the 100 and 010+001 Face 
of Goethite in 0.1'M NaNO, Using a Stern Layer 
Capacitance of 2.23 F/m? (BS Approximation). The site 
density is set at 8 nm? for singly coordinated and doubly 
coordinated surface groups. The calculated overall 6,(pH) 
curve is expressed in mC/m’ BET surface area. The ratio 
of the specific surface area of the two types of faces is 
about 1. 


In Figure 3, the predicted surface charge on the 100 and the 010 
face is given as a function of pH together with the calculated overall 
Particle charge. Note that the pzc of the particle is predicted to be 
around pH 10.0. Another interesting aspect of the predicted charging 
behavior is that, in the region between pH 7.5 and 10.7, the two 
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dominant crystal faces are oppositely charged! The difference in 
charging behavior of the two planes is expected to have an influence on 
the colloid stability of these particles and also on the specific adsorption 
of various ions on these particles. 


Goethite 
a 0.005 M 


+ 0.01M 
° 0.1M 


Surface charge (mC/m? BET) 


FIGURE 4. Predicted Overall 6,(pH) Curves of Goethite (mC/m BET) 
at Three Different NaNO, Concentrations (Full Lines), 
Taking the Weighted Sum of the Gọ-pH Curves of the 
Individual Planes Given in Figure 3. The experimental 
information has been indicated with markers. 


In Figure 4, the predicted charging behavior for goethite at three 
salt levels is compared with experimental data of Hiemstra et al. [22]. 
For each plane, a basic Stern type double layer has been assumed with 
a capacitance of 2.2 F/m’. The Stern layer capacitance was treated as 
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the only adjustable parameter. A very good description of the 
experimental data is obtained using the proton affinity constants as 
predicted by the MUSIC model. Note that the overall pzc is slightly 
dependent on the salt level. This effect is caused by the presence of 
different faces, each with its own double layer, that respond differently 
on a change in salt level. The goethite particle is an example of so- 
called patchwise heterogeneity. The effect of patchwise heterogeneity 
coupled with variable charge surfaces is further discussed by Bolt and 
Van Riemsdijk [36], Van Riemsdijk et al. [37] and Koopal and Van 
Riemsdijk [38]. The high pzc reported here for goethite is in accord 
with other recent findings [39-41]. The pzc reported in literature for 
iron oxides is often in the range from 7 to 10. A pzc lower than 10 for 
a certain iron oxide is not necessarily in conflict with the proton affinity 
constants as predicted by the MUSIC model, because the particle pzc 
depends on the faces being developed and the ratio of the reactive 
groups. Also contamination by CO, may cause the experimental pzc to 
be too low. It is, therefore, of interest to analyze the effect of the 
presence of CO). 


6. ADSORPTION OF CARBONATE SPECIES AND ITS 
EFFECT ON THE CHARGING OF IRON OXIDES 


Dissolved carbon dioxide, carbonic acid and/or (bi)carbonate are 
generally present in natural systems. In laboratory experiments one 
usually tries to exclude carbonate from the system because of the 
possible interference with the measurements. However, for a transfer 
of the laboratory results to the field, it is necessary to know the effect 
of the presence of CO, Only Zachara et al. [18] have reported 
experimentally determined adsorption of carbonate species on 
amorphous iron hydroxide using labelled '*C, but satisfactory modeling 
Was not obtained. The detailed surface structure of amorphous iron 
hydroxide is unknown. Amorphous Fe(OH), can be regarded as water 
rich iron oxide in which the number of oxygens coordinated with iron 
(n) is lower than that of, for instance, hematite (n = 4) and probably 
goethite (n = 3). The coordination number of the surface groups will 
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generally be even lower due to the presence of broken bonds. This 
indicates that the surface groups will be dominated by doubly (n = 2) 
and singly (n = 1) coordinated surface hydroxyls, Fe,OH and FeOH(H) 
respectively. It is assumed here that the surface structure of amorphous 
Fe hydroxide is comparable with that of the reactive faces of goethite 
(010 and 001) which only possess singly and doubly coordinated 
groups, each with a site density of about 8 nm? and that the Stern layer 
capacitance equals that of goethite (C = 2.23 F/m’). 

The relevant protonation reaction for the description of the 
surface charge [22] is given by: 


Fe-OH'® + H, = Fe-0H,™; log K,, = 10.7 (7) 


A species of carbon dioxide is thought to react specifically with 
singly coordinated surface groups [42]. The adsorption reaction can be 
represented as: 


Fe-OH'® + H,CO,’ = Fe-O-“’-C-OOH™ + HO (gy 


log Keii 


or equivalently as: 


Fe-0H"* + HCO, = Fe-0~™-—-C-00H™ + HO (9) 
log K,, = -4.35 i 


Carb 
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0.05 


Total Carbonate = 4.6 umol/| 
[a o 46.4n°/I 
=m 464 m7 


CO, adsorbed (umol/m?) 


pH 
FIGURE 5. Adsorption of CO, on Amorphous Iron Hydroxide for Two 
Different Suspension Concentrations at a Total Carbonate 
Concentration in the System of 4.6 zmol/L. Model 
calculations: solid lines; for parameters see text. Data are 
taken from Zachara et al. [18]. 


It is assumed that the positive charge of the central ions can be 
attributed to one or two of the planes distinguished in the Basic Stern 
model. The charge distribution can be described using the reduced 
valence principle of Pauling [31]. The carbon centre (z = 4+) 
distributes the positive charge over three surrounding ligands, attributing 
to each ligand on the average a valence v of +4/3. The Fe center 
| attributes 1/2+ to the surrounding ligands. This implies that the overall 
‘charge of the surface oxygen will be (1/2 + 4/3 -2 = -1/6) while the 
‘Combined charge of the oxygen and hydroxyl situated at the second 
P ane of the BS model equals -2/6. 
aa The value of the intrinsic log K has been determined using data 
Of Zachara et al. [18]. The data could be described rather well with log 
th = 4.4 (Figure 5). A more traditional formulation of the adsorption 
feaction would be: 
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Fe-OH,”* + HCO; = Fe-OH;,'"’-HCO, (10) 


or 


Fe-0H,™ + HCO, = Fe'*-HCO, + H,O (11) 


With this formulation it is not possible to get a satisfactory 
description of the data. 


CO, adsorbed (umol/n*) 


FIGURE 6. Calculated pH Dependency of the Adsorption of CU, un 
Amorphous Iron (Hydr)oxide for Two Electrolyte Levels 
(0.001 and 0.1 M of a 1:1 Electrolyte) at a Constant CO, 
Pressure. Solid line: -log Poo = 3.5 and dotted line: 
-log Poo: = 2.5. 
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Knowing the affinity of carbon dioxide for singly coordinated 
surface species of iron (hydr)oxides, the expected level of CO, 
adsorption on the surface of iron (hydr)oxides in natural systems can be 
calculated, for example, as function of pH and salt level at various 
levels of a constant P.o, (Figure 6). It is evident that surfaces in 
contact with the atmosphere already may adsorb relatively high amounts 
of carbonic acid, especially at low electrolyte levels. Such a calculation 
indicates that one should be aware of the possibility of a high 
adsorption of carbon dioxide in suspensions used in the laboratory. 


Goethite 


CIP (0.01/0.1 M) 


Total Carbonate (mmol/l) 


FIGURE 7. The Influence of the Presence of Various Levels of Total 
Carbonate in the System on Common Intersection Point 
(CIP) of the Apparent Charging Curves of Goethite 
Suspensions (for 0.01 and 0.1 M) at Two Different 
Suspension Densities. 
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As noted before, the experimentally determined apparent pzc of 
goethites may differ from the theoretical ones. Lower values can be 
due to the influence of CO,. Evans et al. [39] and Zeltner and 
Anderson [41] reported that the careful removal of CO, by purging 
suspensions for quite a long time (2-4 months) increased their common 
intersection point (CIP) by about one pH unit. The influence of CO, 
on the CIP of goethite can now be evaluated. Calculations show that 
only adsorption at the 001 and 010 faces occur. This is in conflict with 
the suggestion of adsorption in the "grooves" at the 100 face as 
proposed by Russell et al. [43] and Zeltner and Anderson [41]. The 
CIP of the apparent charging curves of goethite suspensions at the 0.01 
and 0.1 M electrolyte level have been calculated for various levels of 
total carbonate in the system (Figure 7). 

The calculations indicate that very small amounts of carbonate 
strongly influence the CIP. A goethite suspension stored at pH 8.5 in 
contact with the atmosphere will have a carbonate level of about 3 
mmol/l. Titration of this suspension will give a CIP of about 8.3-8.7, 
depending on the suspension concentration, whereas the theoretical CIP 
without any adsorption of CO, will give a CIP of about 10.1. This 
indicates that the CO, should be removed by purging the suspension. 
However complete removal of CO, by purging is a very slow process 
if a suspension is only purged without any acidification, as is frequently 
done. This is illustrated in Figures 8 and 9. A volume of 50 ml of a 
10° and 10* M NaHCO, solution in a closed vessel (100 ml) was 
purged with purified moist N, gas (100 ml/min). 

Figure 8 shows that the pH rises very quickly during the first 
hours of the experiment, and Figure 9 shows that the corresponding 
removal of CO, is rather fast. However, increasing pH decreases the 
Peo Strongly, causing a strong reduction of the flux of removed CO). 
It can be concluded that with such a purging process even after two 
days less than 50 % can be removed. The resulting CO, levels may 
still reduce the CIP by about 0.5-1 pH unit. The removal of CO, can 
be enhanced by keeping the pH constant during the process. From a 
10° M NaHCO, solution (pH 8.1) all CO, can be removed in this 
manner in approximately 24 hours (Figure 9). Figure 7 indicates that 


FIGURE 8. 
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| 0.001 M NaHCO, 


25°C 
50 ml solution 
100 mi NJmin 
o 20 40 60 
time (hours) 


Experimentally Determined pH Development During the 
Removal of CO, from a 10* and 10° M NaHCO, (in 0.05 
M NaNO,) Solution by Purging It with Moist Purified N, 


even a very small amount of carbon dioxide will change the CIP 
Strongly. It implies also that all solutions added to purged suspensions 
Should be free of CO,, which is quite difficult for NaOH solutions. 
The shift of the CIP of goethite from 10.1 to lower values is 
‘Mainly due to the titration of free carbonate species in solution. The 
Adsorption of CO, as such hardly influences the surface charge. 
‘Calculations show that the electric potential in the surface layer and the 
‘Second adsorption layer is hardly influenced by the adsorption. It 
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0.001 M NaHCO, 
pH constant 


0.0001 M NaHCO, 
pH free 


Percentage CO, removed 
8 


time (hours) 
FIGURE 9. The Removal of CO, from a N, Purged 50ml 10* and 10° 
M NaHCO, (in 0.05 M NaNO,) Solution Without pH 
Control (pH Free). Calculated from the pH change given 
in Figure 8 and from 50 ml of 10° M NaHCO, (in 0.05 M 
NaNO,) solution in which the pH was kept constant at 8.15 
by titrating acid. 


implies that the isoelectric point (IEP) of the goethite is not changed 
upon adsorption of moderate amounts of carbonic acid. This prediction 
is in accord with the experiments of Zeltner and Anderson [41] who 
found a CO, independent IEP whereas their CIP in the titration curves 
was strongly CO, dependent. It can be concluded that the pzc of metal 
oxides with a high pzc like goethite can only be accurately established 
with very careful experimental work. 

In several cases the reported values of the pzc of goethites are 
below pH 8.0. This cannot be completely understood from the above 
given analysis of the influence of CO, on titration curves. Other factors 
which may influence the CIP are the presence of silica [44] or other 
impurities at the surface [45] and also a variation in the relative 
presence of crystal planes or surface composition. 
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7. ADSORPTION MODELS FOR SOIL ORGANIC 
MATTER 


Humified organic matter, e.g. humic acids, fulvic acids, and 
humin, in the natural environment may be a very important sorbent for 
metals in soils, sediments and aquifers. These materials are 
characterized by a pH dependent negative surface charge and surface 
potential. The surface charge is caused by dissociation of various types 
of functional groups. There is a wide spread in intrinsic proton affinity 
constants for the different functional groups, ranging from around log 
K = 2 to log K = 11. The interaction of ions with humified organic 
matter is influenced by the chemical heterogeneity and by the variable 
surface potential. The variable charge character of these materials 
shows up clearly in the salt level dependency of the charging curves. 
The charging curves are equivalent with plots of the degree of proton 
dissociation as a function of pH. For a proper physical chemical ion 
adsorption model, it is necessary to account separately for the effect of 
the chemical heterogeneity and the effect of the variable surface 
potential on ion adsorption. For the assessment of the surface potential 
a double layer model is required. It is thus necessary to establish an 
appropriate double layer model and assess the chemical heterogeneity 
before one can construct a realistic adsorption model. All models that 
have been put forward so far make several a priori assumptions both 
with respect to the electrostatic part and with respect to the treatment 
of the chemical heterogeneity. Most often the electrostatic part is not 
even dealt with explicitly, i.e., that the electrostatic effects are 
incorporated in the chemical heterogeneity. Recently, De Wit et al. 
[46,47] have proposed a procedure that allows for the assessment of the 
double layer model and its parameters independently from the 
estimation of the intrinsic chemical affinity distribution. The principle 
Of the method is that it is possible to filter the electrostatic effect out 
Of the data by plotting the charging curve not as a function of pH but 
aS a function of pH,. It follows from Equation (3) that H, is equal to 
Me proton concentration in the bulk of the solution times a Boltzmann 

actor. The Boltzmann factor can be calculated using a double layer 
aodel. When a proper double layer model is used, the charging curves 
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at different ionic strength values should all coincide into one curve 
when plotted as a function of pH,. We have called this curve the 
master curve. From this master curve the intrinsic chemical 
heterogeneity, or the intrinsic affinity distribution, can be estimated 
using special techniques [48-51]. A serious problem in the affinity 
analysis is the treatment of experimental error in the data. Advanced 
methods are necessary to obtain sensible results [52]. It should be 
realized that information is lost when experimental error is present. 
When the data handling is inappropriate many spurious peaks will 
appear in the affinity spectrum [48,52], whereas, with a proper data 
treatment, the fine structure of the affinity spectrum cannot be 
recovered due to the uncertainty in the data. Applying the master curve 
concept to charging curves of humic and fulvic acids [47] shows that 
a spherical double layer model is in some cases appropriate, whereas in 
other cases, a simple Gouy Chapman double layer model suffices. For 
small spherical entities the salt effect on the charging curves is less than 
for larger entities. If the radius is larger than ca. 10 nm the spherical 
model becomes equivalent to the Gouy Chapman model. 


8. APPLICATION TO NATURAL SYSTEMS 


Doth the model development and the collection of experimental 
at present, not far enough advanced to allow for a direct 
in of the behavior of metals in natural systems based on 
chemical models only. The knowledge with respect to 
mn phenomena on metal (hydr)oxides is far more advanced than 
for humic materials. 

he variability in the extent of sorption for different soils, 
s, or aquifer materials originates from differences in~solid 
aracteristics and from differences in parameters like pH and 
tential. The physical chemical models allow the construction 
mework that allows for the "normalization" of the wide 
y that is observed in nature. This is a very important aspect 
‘velopment of quality criteria for soil and sediment [53,54]. 
ibilities for transfer of knowledge from relatively pure model 
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systems to natural systems will be illustrated using orthophosphate as 
an example. Phosphate has the advantage that there is no direct 
interaction between this ion and humic materials. In many natural 
systems the iron and aluminum (hydr)oxides are the dominant reactive 
surfaces for phosphate sorption. Another big difference between metal 
ion adsorption and phosphate adsorption is that the effect of pH on 
phosphate adsorption is relatively small. 

From adsorption studies with pure crystalline metal (hydr)oxides, 
e.g. goethite, it is observed that the adsorption isotherm is very flat 
above a few tenths of a millimole phosphate [55]. With gibbsite as 
sorbent it is found that the sorption reaction may proceed far beyond 
the adsorption maximum [56]. This was interpreted as the formation 
of coatings of metal phosphate at the expense of the underlying metal 
hydroxide. Based on these ideas, a kinetic diffusion-precipitation model 
was developed that is applicable for an arbitrary mixture of reactive 
particles that may differ in size, shape and chemical composition 
[57,58]. The total phosphate sorption in soil is thus a combination of 
reversible adsorption and diffusion-precipitation. The suggestion that 
phosphate sorption in soil may occur through the formation of metal 
phosphate coatings on metal (hydr)oxides has recently been confirmed 
for a wide range of soils using surface analysis techniques [59]. For 
both processes the amount of iron plus aluminum extractable with 
ammonium oxalate has been used to estimate the (ad)sorption capacity 
[60,61]. The reversible part of the total sorption in a given soil sample 
can be estimated using a newly developed desorption technique [62]. 
It is interesting to note that the shape of the reversible adsorption 
isotherms for the soils considered is quite similar to the shape of the 
phosphate adsorption isotherm for metal (hydr)oxides as discussed 
above. In contrast, the shape of the total sorption versus the phosphate 
concentration shows a considerable slope above a few tenths of a 
‘millimole. Transport of phosphate through soil-filled columns could be 
satisfactorily predicted [63] using a reversible adsorption model in 
Combination with the diffusion-precipitation model with independently 
estimated model parameters. 
The example discussed above illustrates that the physical 
[chemical frame-work that originated from work with model systems 
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was extremely powerful as a basis for the development of models that 
describe the behavior of phosphate in soils. It also illustrates that, apart 
from reversible adsorption on reactive surfaces, other processes may be 
of importance in determining the fate of chemicals in the environment. 


9. SUMMARY 


Adsorption of ions in natural systems is quite complex and at 
present not completely understood. The complexity is due to the fact 
that natural systems are mixtures of various reactive surfaces. This in 
itself leads to chemical heterogeneity that affects ion adsorption. An 
extra complicating effect is that the sorbent components present in 
nature, like metal oxides or natural organic matter are heterogenous in 
itself. 

The combination of various types of heterogeneity, with variable 
charge characteristics and competition between various ions for the 
same sorption sites are the main ingredients for the complexity of the 
situation. Some mechanistic ion adsorption models for 
metal(hydr)oxides and natural organic matter are discussed. It is shown 
that carbonate contamination may easily lead to a too low value of the 
experimentally determined pzc of metal(hydr)oxides with a high pzc. 
It is shown why this error may still be significant even when the system 
is flushed with nitrogen gas. In conclusion, it is discussed how insights 
derived from fundamental studies can be applied to natural systems. 
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